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Geometrical spin manipulation in Dirac flakes 
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We investigate numerically the spin properties of electrons in flakes made of materials described 
by the Dirac equation, at the presence of intrinsic spin-orbit-coupling(SOC). We show that electrons 
flowing along the borders of flakes via edge states, become helically spin-polarized for strong SOC, 
for materials with and without a gap at the Fermi energy, corresponding to the massive and massless 
Dirac equation respectively. The helically spin-polarized electrons cause geometrical spin splitting 
on opposite sides of the flakes, leading to spin-resolved transport controlled by the flake’s geometry 
in a multi-terminal device setup. A simple analytical model containing the basic ingredients of the 
problem is introduced to get an insight of the helical mechanism, along with our numerical results 
which are based on an effective tight-binding model. 

PACS numbers: 72.80.Vp, 71.70.Ej, 75.76.+J, 73.63.Kv 


INTRODUCTION 


Mesoscopic low-dimensional materials described by the 
Dirac equation for relativistic particles is a rapid de¬ 
veloping field providing a vast area of fundamental re¬ 
search and promising a wide range of applications. One 
of the earliest examples is graphene, a two-dimensional 
(2D) sheet of carbon atoms arranged in a honeycomb lat¬ 
tice structure resulting in an electronic behavior at the 
Fermi level that resembles relativistic massless particles 
described by the massless Dirac equation [l:, 2}. Another 
more recent example is transition metal dichalcogenide 
(TMD) monolayers , with honeycomb lattice struc¬ 
ture similar to graphene but with a varying gap at the 
Fermi level, for instance 1.35 eV for WS 2 and 1.83 eV 
for M 0 S 2 . Alternatively a considerable gap of up to 
0.278 eV can be created in graphene by embedding it 
on boron nitride (BN) substrates [l^ - [l3 | . Graphene on 
BN substrates and TMD monolayers can be effectively 
described by the massive Dirac equation, instead of the 
massless one which describes pristine graphene. 

Dirac materials can be easily formed in confined dot¬ 
like structures of various sizes and shapes, known as 
flakes, having been studied theoretically and demon¬ 
strated experimentally jlj, Il5l - f26| . Specifically for TMD 
flakes the trigonal shape is a natural state in their exper¬ 
imental preparation [2 51I2(jf. The confinement in Dirac 
flakes leads to unconventional electronic properties such 
as the edge states, whose electronic density is concen¬ 
trated at the edges of the flakes 16 jlil |2 t 1 leading 
to novel topological phenomena 28 - 3lj. Edge states 
are fundamentally a consequence of the honeycomb lat¬ 
tice morphology present in Dirac materials that favors 
destructive interference effects when zig-zag edges are 
present. 

On the other hand spin-orbit-coupling(SOC) in Dirac 
materials offers the possibility to manipulate the elec¬ 
tron spin for integration in spintronic applications 32 - 
[3^ |. Despite the weak SOC in graphene, the interaction 
with the system’s confinement, particularly in graphene 


nanoribbons, has been shown to produce novel quantum 
phases like the Quantum-Spin-Hall(QSH) effect 36-381]. 
The QSH effect has sparked a continuously increasing in¬ 
terest in exotic materials that display related behaviors 
known as topological insulators. Also spin-transport ef¬ 
fects have been investigated in graphene flakes via DFT 
simulation[ 3 ^, [Icj. In the same sense, TMD monolay¬ 
ers and other emergent 2D materials are also promis¬ 
ing candidates for relevant theoretical investigations fll] 
and applications since the SOC strength can be much 
higher than graphene, leading to large spin-spilittings up 
to 0.46 eV in some cases of TMD monolayers [42| . 


In the current work we analyze the spin properties 
of electrons in Dirac flakes with intrinsic SOC via the 
numerical calculation of the transmission probabilities 
along the border of the flakes. We show that for suffi¬ 
ciently strong SOC the electrons that propagate in op¬ 
posite directions along the flake’s border via the edge 
states, have opposite spins, leading to helically spin- 
polarized electrons. We demonstrate this property in a 
multi-terminal device setup naturally formed by attach¬ 
ing three leads at the corners of a trigonal flake. Ap¬ 
propriate tuning of the chemical potential of the leads 
creates spin-resolved electronic transport with opposite 
spin-polarization along the opposite sides of the trian¬ 
gle. This allows a geometrical manipulation of the spin 
intrinsically without the need to apply external fields. 
The helically spin-polarized electrons are present in ma¬ 
terials with and without a gap at the Fermi energy, de¬ 
scribed by the massive and massless Dirac equation, re¬ 
spectively, which could correspond to TMD monolayers 
and graphene. Additionally, we introduce a simplified 
analytical model, that contains the fundamental ingredi¬ 
ents of the numerical model, allowing us to get a better 
understanding of our numerical results. 
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FIG. 1. A trigonal Dirac flake with the characteristic honey¬ 
comb lattice structure consisting of two sublattices A and B. 
An example of the spin-orbit coupling interaction involving 
spin-dependent hopping amplitude between second-nearest 
neighbor sites can be seen inside a single hexagon. Three 
leads can be attached naturally to the corners of the trigonal 
flake. The arrows across the flake’s border indicate the he¬ 
lically spin-polarized electronic flow, satisfying the relations 
for the transmission probabilities between different terminals, 
shown on the right. We characterize the flake size by the 
length L of the triangle base, in units of the lattice constant 
a. The edge state mechanism due to the honeycomb lattice 
structure can be seen on the left along with a schematic show¬ 
ing how to exchange the A and B atoms by cutting the flake 
in different orientations. 


MODEL 


In this section, we present our numerical tight-binding 
model to simulate flakes made out of Dirac materials at 
the presence of SOC. Dirac materials can be described 
by the massive Dirac Hamiltonian 41, 43l-47| 


Ho = hvf{k x T z ai + k y a 2 ) + Va 3 


(1) 


where Vf is the Fermi velocity of electrons whose value 
depends on the material under investigation, and cq be¬ 
ing the Pauli matrices acting on orbital space. Sym¬ 
bol t z = ±1 denotes the non-equivalent valleys that 
are present in graphene and other Dirac materials like 
TMD monolayers, at the six corners of their hexagonal 
Brillouin zone. Eq. |T] describes relativistic particles of 
mass V, while the speed of light c is replaced by Vf. 
Different values of V classify different materials for ex¬ 
ample V = 0 meV and Vf ~ 10 6 m/s corresponds to 
graphene, while finite V could describe TMD monolay¬ 
ers and graphene on BN substrates 14. 41. 43l . 44 1. 

In order to perform our numerical calculations, we con¬ 
sider an effective tight-binding model, consisting of a hon¬ 
eycomb lattice with first nearest neighbor hopping, along 
with an onsite staggered potential simulating the mass 


term in Eq. [TJ 

Ho = ^2 e i c l A t c i/j + ^2 + c l/i c b) ’ (2) 

i (*b> 


where c \ (cj^) is the creation (annihilation) operator for 
spin /r at site i while = 1 eV is a uniform hopping 
between all the nearest neighbor lattice sites. The on¬ 
site potential e ? ; is V and —V on A and B sublattice 
sites respectively with V = 0 corresponding to graphene. 
The staggered potential breaks the inversion symmetry 
resulting in a gap 2V at the Fermi energy for infinite 
unbounded systems. At low energies the effective tight- 
binding Hamiltonian Eq. [2] transforms to Eq. [l] with 
hvf = t '^ a where a is the lattice constant. This effec¬ 
tive tight-binding model can be thought as a numerical 
version of the massive Dirac equation. 

The SOC can be introduced in the tight-binding model 
by considering spin-dependent hoppings a long ne xt near¬ 
est neighbors in the honeycomb lattice as[36l 1431] . 

H S oc=a ^2 ( 3 ) 

where the sum runs over next nearest neighbors on 
all the lattice sites. The spin-dependent amplitude is 
Vij = 1(—1) when the electron at site i makes a anti- 
clockwise(clockwise) turn in order to hop from the first 
nearest to the second nearest neighbor j. represented by 
the curved arrows inside the honeycomb lattice of the 
trigonal flake in Fig. [1] The strength of the intrinsic 
SOC is determined by a. We note that the SOC inter¬ 
action described by Eq. [3] involves connections between 
sites belonging to the same sublattices A or B, in con¬ 
junction with the first nearest neighbor term in Eq. [2] 
which involves connections only between A and B sub¬ 
lattice sites. Eq. [3] is responsible for the QSH effect in 
confined graphene systems, acting as a bridge between 
graphene and topological insulators. The total Hamilto¬ 
nian of our system is 

H = H 0 + H SOC ■ (4) 

By applying hard-wall boundary conditions on Eq. [4] we 
can simulate the trigonal flakes, like the one shown in 
Fig. [Q 

In order to calculate the transmission probabilities via 
the edges of the trigonal flake, we attach perfect semi¬ 
infinite linear chains described by Hid = X^( c l/i c i+1 m + 
c l+ w ith hopping t=leV at it’s three corners as 

shown in Fig. [I] forming this way a three terminal de¬ 
vice. The energy E, the SOC strength a and V, are all 
reported in units of t. We observe that Eq. [3] contains 
only s z so that the total Hamiltonian of the system Eq. 
[|] can be split in spin up and down diagonal blocks, re¬ 
sulting in zero transmission probability between opposite 
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FIG. 2. a) The energy levels of a trangular Dirac flake for 
L=8, V=0.4 and different SOC strengths a near the gap in¬ 
duced by V. Fourteen edge states at E=-V can be distin¬ 
guished for a = 0 due to sublattice symmetry, dispersing 
inside the gap for finite a. The levels become almost homoge¬ 
neously distributed for large a. b) The respective wavefunc- 
tion for a = 0 at E=-V is an edge state being non-zero only 
on B sublattice sites, c) The edge state is enchanched by the 
SOC becoming completely concentrated at the edges, d) The 
wavefunction is accumulated at the corners for large a. 

spins(see Appendix). Therefore the transmission proba¬ 
bility between different terminals can be written as Tj'- 
where i,j=l,2,3 runs over the terminals and /i 4- de¬ 
notes the spin orientation of the electron. We define also 
the spin transmission as = tJ- — T-j which can be 
used to estimate the degree of spin-polarization. TT can 
be calculated by using the Green’s function £t, E) = 

{E+ -H- T, 1 {E+) - E 2 (E+) - T, 3 (E + ))~ 1 with E+ = 
E + ir /, where E is the incident energy and r) —> 0 + , while 
E i(E) is the self energy of each semi-infinite chain used 
as leads, being diagonal in the spin-basis, since we ne¬ 
glect the SOC inside the leads(see Appendix). By using 
the velocities inside the leads according to the Fisher-Lee 
relations we get Tfi = (4 - E 2 )\G(j, i, //, E)\ 2 . 

EDGE STATES 

In this section, we analyze the edge states present 
in the trigonal zig-zag flakes made of Dirac materials, 
which are responsible for the electronic behavior near the 
Fermi level. In the lattice representation, the number of 
edge states can be derived by extending chiral symme¬ 
try arguments applied on the honeycomb lattice which is 
bipartite[48(. In general the Dirac flakes can be split into 
two sublattices with Na number of A and Nb number 
of B sublattice sites, as shown in Fig. [L] The type of 
the outermost atoms at the zig-zag edges of the flake de¬ 
termines which number, either Na or Nb is larger. For 


the configuration in Fig. [l] with B type atoms, we have 
Nb > Na , while Nb — Na = L — 1, with L being the 
length of the triangle basis in units of the lattice con¬ 
stant a. We note that the type of atoms A and B can 
be easily exchanged by changing the orientation of the 
triangle base on an infinite graphene sheet by 60 degrees, 
compared to the base of the flake in Fig. [TJ resulting in 
Na > Nb in this case. Since the spin-orbit coupling given 
by Eq. [3] contains only s z it does not mix the spin up 
with the spin down orientations, so that the total Hamil¬ 
tonian of the system can always be split into spin up and 
spin down diagonal blocks. Consequently we can write 
down the Schrodinger difference equations centered on A 
and B atoms for spin up (down) denoted by /i = 1(—1) 
as, 

(■ E ~ e A )^A,i = * Ej *B,j + OL Vij^Aj 

(5) 

{E - e B )K,i = * E j Kj + ^ Ei B,j 

where a(b),i is the wavefunction amplitude on site i 
of each sublattice and (-a(b) is the staggered on-site po¬ 
tential. In both sets of equations, the first term on the 
right-hand side comes from the nearest neighbour hop¬ 
ping, while the second term is associated with the second 
nearest neighbor hopping with spin dependent amplitude 
Vi.j. as in Eq. [3] In the absence of SOC (a = 0), we can 
easily derive the number of edge states by algebraic ar¬ 
guments applied on Eq. 0 For e A = V,6 b = —V the 
equations become, 

(E-V^^t^Kj ( 6 - 1 ) 

3 

(E+V)^ B .=tY,n,r ( 6 - 2 ) 

3 

Eq. 16.11 is a set of N A equations with Nb unknowns 
written on the A sites, while Eq. I6.2l is a set of Nb equa¬ 
tions with Na unknowns written on the B sites, respec¬ 
tively. We observe that for E=-V, Eq. 16.21 transforms to 
the homogeneous set 0 = t Efc x ^a k which can only be 
satisfied by assuming that 4/^ = 0 since there are more 
equations than unknowns. On the other hand Eq. 16.11 
transforms to -2WF^ = tJ2 k ^B,k =► 0 = *E k^B,k 
which gives but results also in Nb — Na linearly in¬ 
dependent solutions, since there are less equations than 
unknowns. 

Consequently for a = 0 there are at least Nb — N a 
states at E=-V with non-zero amplitudes only on the B 
sublattice sites. For V=0 we derive the well-known case 
of bipartite lattices where there are at least Nb — N A 
states at E=0 due to the chiral symmetry . With the 
application of the staggered potential (V ^ 0) the energy 
of these states is shifted to E=-V. Considering the spin 
these states become doubly degenerate. 
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Due to the destructive interference effects favored by 
the honeycomb lattice when zig-zag edges are present, 
these states either at E=-V or E=0 are edge states with 
their corresponding wavefunctions concentrated at zig¬ 
zag edges of the flakes. The edge state mechanism can 
be robustly understood by isolating part of the zig-zag 
edge as shown in Fig. |T] at the right side of the flake 
schematic, where there are three B sites surrounding one 
A site. The edge states can be visualized as two incoming 
electronic waves from two B sites 1 and 2, interfering and 
giving an outgoing wave at site 3 site via the A site. For 
E=0 there is destructive interference of the two incoming 
waves due to a phase difference between the wavefunction 
amplitudes on the respective sites 1 and 2 , introduced 
by the Bloch’s theorem, as can be seen easily in a semi¬ 
infinite graphene sheet with one zig-zag edge 27]. These 
two waves from sites 1 and 2 interefere destructively via 
the central A site whose wavefunction is zero, giving zero 
amplitude at site 3 for states at the ends of the Brillouin 
zone where the phase difference is 7r or a reduced am¬ 
plitude otherwise. This intereference process continues 
inside the rest of the honeycomb lattice and is present 
whenever zig-zag edges exist at the borders of a confined 
graphene structure. The application of the staggered po¬ 
tential, does not alter this basic topology of the honey¬ 
comb lattice structure which is responsible for the edge 
states, since it gives a potential -V on all the B sites 1,2,3 
and potential V on site A, resulting in the shifting of the 
edge state mechanism to energy E=-V. We note that the 
intereference mechanism responsible for the edge states 
vanishes for the armchair edge morphology. 

According to the above algebraic arguments, the en¬ 
ergy of the edge states can be shifted at the opposite 
(conduction) side of the energy spectrum at E=V by 
simply exchanging the potential on A and B sublattice 
sites, so that the outermost sites at the zig-zag edges in 
Fig. |T| have on-site potential V instead of -V. We sum¬ 
marize that for Dirac flakes with zig-zag edges, in the 
absence of SOC, there are at least | Na — Nb | edge 
states at E=V(E=-V) when the potential of the outer¬ 
most atoms of the zig-zag edges is V(-V). Therefore, there 
are always edge states whose energy is determined by 
the type of atoms at the zig-zag edges of the flake, while 
their number increases linearly with the system size L as 

Na — Nb \= L — 1 . In realistic systems, our results 
imply that the energy of the edge states can be tuned 
near either the valence or the conduction band edge, by 
appropriately cutting the Dirac flakes, in such a way that 
the type of atoms A and B are exchanged. 

When the SOC is introduced in the system, for a ^ 0 
in Eq. [5j the edge states disperse inside the gap created 
either, due the staggered potential V, or due to the fi¬ 
nite flake size. For instance, this phenomenon has been 
observed in M 0 S 2 flakes by multi-orbital tight-binding 
simulation [24|. In our case this behavior can be under¬ 
stood by the SOC term in Eq. [5] which alters the values 


of E for which these equations can be satisfied accord¬ 
ing to the algebraic arguments we presented. The energy 
E=-V is no longer a solution of these equations, while the 
wavefunction amplitude becomes gradually finite on both 
sublattice sites A and B as a is increased. In the simu¬ 
lation of specific materials the SOC term Eq. [3] has to 
be applied only on one of the sublattices. If it is applied 
only on A then the results we derived for a = 0 remain 
valid, since the wavefunction amplitude on A sites is zero 
and any perturbation applied on this sublattice will not 
affect the system’s properties. 

We have verified our analytical results numerically. An 
example can be seen in Fig. [2^ where the energy levels 
for a trigonal flake with L=8 and V=0.4 are shown, near 
the energy gap. All states come in pairs due to the spin 
degeneracy, which is preserved even for for finite a. In 
agreement with our algebraic arguments there are 2 (Nb — 
Na) = 14 edge states at E=-V, which disperse inside the 
gap for finite a. For sufficiently large a the states tend 
to repulse, becoming almost homogeneously distributed 
inside the whole gap, despite being initially concentrated 
at it’s lower end(E=-V). In general a significant interplay 
between the SOC and the gap generated by the staggered 
potential is expected. Additionally the edge states due 
to zig-zag edges become gradually mixed with the edge 
states created due to the SOC in analogy to the edge state 
mechanism in the Quantum-Hall effect. In this sense, the 
edge state mechanism originating from the honeycomb 
lattice structure is enhanced by the SOC, however it’s 
contribution becomes irrelevant for strong SOC. 

Some examples of these edge states can be seen in Fig. 
Eb,c,d where the wavefunction propability is shown, rep¬ 
resented by the radius of circles at each site. For a = 0 
the amplitude of the edge state at E=-0.4, is non-zero 
only on the B sites, with the largest amplitude residing 
along the edges. There is almost perfect concentration 
of the wavefunction at the zig-zag edges of the flakes for 
finite a. A gradual accumulation at the corners of the 
triangle can be observed for large SOC (a = 0.5). This 
is due to the different lattice connectivity between the 
sites at the corners of the triangle and the ones along 
the zig-zag edge. Particularly the corner sites are con¬ 
nected via two spin-dependent hoppings to second near¬ 
est neighbors, instead of four for sites along the zig-zag 
edge, coming from Eq. [3] Also, each of the corners can 
be thought as the intersection between two zig-zag chains 
causing additional intereference effects, that cannot be 
distinguished for a = 0 since the wavefunction probabil¬ 
ity is zero on A sites. We note that the results are similar 
for the spin up and down cases. 


TRANSMISSION PROBABILITIES 

In the current section we investigate the transmission 
probabilities via the edges of the trigonal Dirac flakes 
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V=0 


V=0.4 


V=0.8 



oc=0.05 


a=0.20 


oc=0.50 



FIG. 5. T 12 for a = 0.5 and different flake sizes L. The helical 
regime is retained for all L independently of V. 


FIG. 3. The transmission probabilities T^ 2 between terminals 
1 and 2 for spin up and down, different a and V. For strong 
SOC a = 0.5 a significant difference between up and down 
is formed indicated by an arrow where T± 2 ~ 1 and T^ 2 « 0, 
signifying helically fully spin-polarized electrons flowing along 
the edge. The helical regime is shifted along the energy for 
increasing V. 



with SOC. In order to probe the transmission we form 
a three terminal device setup by attaching linear chains 
at the three corners of the triangle, as shown in Fig. 
□ In Fig. [3] we show the spin-dependent transmission 
probabilities T^ 2 for L=8 along one side of the triangle 
from terminal 1 to 2 for the spin-up and spin-down cases. 
We have performed the calculation for different values of 
the staggered potential V and the SOC strength a. The 
electronic flow inside the energy area where we plot our 
results, is carried via edge states like the ones in Fig. [2] 
for finite a. For weak SOC (a = 0.05) resonant trans¬ 
mission via the respective energy levels of the system 
can be seen as sharp peaks. Increasing a causes broad¬ 
ening of these peaks resulting in the smooth fluctuations 
shown. Clearly, a significant difference between the oppo¬ 
site spins is starting to show for sufficiently large a. For 
a = 0.5 a small regime is formed for every V, indicated 
with a small arrow, where there is almost perfect propa- 


L i 2 ~ 0). We 


gation of solely spin-up electrons (T{ 2 « 1 , T-j 
note that large SOC interactions have been demonstrated 
in TMD monolayers 42] which correspond to finite V in 
our model. 



E 



izvj _i_i_i_i_i— 

0 0.2 0.4 0.6 


FIG. 4. a) Density color plot of the spin transmission versus 
E and a. Inside the yellow area where T{ 2 * 1 the electrons 
are helically spin polarized. Larger a is required for materials 
with wider gaps(larger V)in order to observe the helicity. b) 
Some respective cases of T{ 2 versus E for different a. 


This behavior can be attributed to two mechanisms. 
One is the helical currents at the edges induced by the 
QSH effect followed by the Laughlin’s argument, in anal¬ 
ogy with the zig-zag nanoribbons with SOC[2fj]. An ad¬ 
ditional mechanism can be revealed by a carefull obser¬ 
vation of the zig-zag edge. In Fig. |T]we can see that the 
electronic flow from terminal 1 to 2 is carried through 
the outermost B (black) sites via anticlockwise turns cor¬ 
responding to positive spin-flip amplitude v = +1 in Eq. 
[31 implying spin-up polarization. The other possible flow 
via A(red) sites across the edge corresponding to spin- 
down polarization is partially supressed since the wave- 
function amplitude is generally smaller than on B sites, 
being a remnant of the chiral symmetry for a = 0, as we 
have shown in the previous section. 

The spin-polarized regime is shifted slightly by increas¬ 
ing V, due to the scattering via the positive potential V 
encountered by the incident electrons from the leads at¬ 
tached to the corners of the triangle. Potential V can be 
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thought as a barrier that the electrons need to overcome, 
by acquiring sufficiently large incident energy E. We have 
verified that the transmission probabilities between 

different terminals satisfy the conditions , which 

are a consequence of the triangle’s symmetry. Combined 
with Tl 2 « 1, 2 « 0 they lead to the conditions shown 

at the right side of Fig. [I] As long as these conditions 
are satisfied the electrons propagating in opposite direc¬ 
tions along the edges of the flake, will have also opposite 
spins. In this sense, the electrons in the trigonal Dirac 
flakes become helically spin-polarized along the edge for 
sufficiently strong SOC. Additionally, we observe that 
the transmission is enchanched by the SOC due to the 
accumulation of the wavefunction at the corners of the 
triangle where the leads are attached, as can be seen in 
Fig. [2ji. However the helicity supresses the electronic 
flow for one of the spin orientations and consequently 
the total (charge) transmission remains almost one for 
every a at energies where the transmission occurs via 
the energy levels of the flake. 

To further analyze the helical regime, in Fig. [4]we show 
a density color-plot of the spin transmission T® 2 versus 
E and a, along with individual cases for constant a in 
the bottom. The helical(yellow) regime can be clearly 
distinguished where Tf 2 becomes maximum (T-f 2 = 1), 
signifying fully spin-polarized electrons that propagate 
helically along the flake’s edge. We notice that for larger 
V, corresponding to materials with wider gaps, larger 
SOC is needed in order to observe the helically polarized 
electrons. This is due to the interplay between the gap 
generated by _the staggered potential and the one due 
to the SOC[36j. In order for the QSH effect to manifest, 
being essentially responsible for the helicity, the SOC has 
to dominate the gap over the staggered potential. The 
helical regime is retained for all flake sizes, ranging from 
tenth to a few hundred atoms, as can be seen in Fig. El 
with more fluctuations created, due to the denser energy 
levels as the size is increased. 

By appropriately tuning the chemical potentials fj,i 
corresponding to each terminal i in the three-terminal 
setup[49(, spin-polarized currents with opposite polariza¬ 
tion flow along the opposite sides of the triangle. For 
Mi = M3 spin currents If = — If = ^t(mi — M 2 ) flow at ter¬ 
minals 1 and 3 with opposite spin polarizations. There¬ 
fore, spin-resolved transport along the different sides of 
the flake is possible, allowing in this sense a geometrical 
manipulation of the electron spin, since the path of the 
spin-polarized electrons can be controlled by the flake’s 
shape. 

We expect similar phenomena for differently shaped 
flakes such as hexagonal, which could offer additional 
possibilities to control the spin dependent transmission 
by attaching more leads at the flake’s corners. Also im¬ 
purities or weak disorder should not affect our results 
since the system is protected from backscattering due to 



FIG. 6. a) A simple model simulating the three-terminal 
setup of the trigonal Dirac flake, the edge state mechanism 
along with the influence of the intrinsic SOC. b) T( 2 for differ¬ 
ent q and V. The curves describe the trend of the numerical 
results for all cases. 


the time reversal symmetry, as in QSH systems. 

In order to further investigate the helical mechanism 
we introduce a simple model consisting of a triangle with 
three sites with constant on-site potential V in analogy 
with the potential at the corners of the flake in Fig. |T| 
and three linear chains attached to it’s corners acting as 
leads, as shown in Fig. [6ji. Also, we consider spin depen¬ 
dent hoppings of strength a inside the triangle simulating 
the SOC. This model can be considered as an isolation 
of the SOC mechanism for each of the sublattices A or 
B (Eq. [3| inside one hexagon in the honeycomb lattice, 
corresponding to L=1 for the trigonal Dirac flake. Alter¬ 
natively it can be thought as a robust simulation of the 
numerical model, since it contains the basic ingredients 
of the problem, namely the fact that the electrons prop¬ 
agate in a straight orbit along the border of the trigonal 
flake due to the edge states, under the influence of the in¬ 
trinsic SOC. We note that due to the way we oriented the 
triangle the transmission from site 1 to 2 T[ 2 takes nega¬ 
tive values. In Fig. [5 }d we show the spin transmission Tf 2 
for different V and SOC strengths. We show —T( 2 in or¬ 
der to lie in the positive axis. The calculation details are 
presented in the Appendix. The analytical curves catch 
the trend of the numerical results for the trigonal flake 
with maximum T( 2 « 1 for sufficiently large a, which is 
shifted across the energy for finite V. 


CONCLUDING REMARKS 

We have presented a numerical study of the spin prop¬ 
erties of electrons in flakes made of Dirac materials at 
the presence of intrinsic spin-orbit-coupling (SOC). A de¬ 
tailed analysis of the edge states is also presented, based 
on sublattice symmetry arguments. The flakes are sim¬ 
ulated by an effective tight-binding model, which can be 
thought as a numerical version of the Dirac equation. We 
probe the spin dependent transport by forming a multi- 
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terminal device which allows the calculation of the trans¬ 
mission probabilities along the edges of the flakes via the 
Green’s function formalism. 

We have shown that the electrons propagating along 
the border of the flakes via the edge states, become heli¬ 
cally spin polarized for strong SOC, for both the massless 
and massive cases, corresponding to Dirac materials with 
and without a gap at the Fermi energy, respectively. The 
helicity can be easily detected in a three-terminal setup 
formed naturally by attaching Id leads at the corners of 
a trigonal flake. Appropriate tuning of the chemical po¬ 
tential of each lead, creates spin-resolved transport with 
opposite polarization between the opposite sides of the 
flakes, allowing a geometrical manipulation of the elec¬ 
tron spin. 

In other words, we have shown that spin-resolved 
transport can be created intrinsically in flakes made of 
Dirac materials by forming multi-terminal devices with¬ 
out the need of external fields, the only requirement being 
sufficiently strong intrinsic SOC. We hope that our work 
will motivate further investigation of the SOC effects in 
2 d materials in conjunction with topological effects due 
to confinement in flakes and other nano-structures. 
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Appendix: Simple model 

In the following Appendix, we present the derivation of 
the transmission probabilities for a simple system consist¬ 
ing of a triangle with three sites, and three linear chains 
attached to it’s corners acting as leads. We consider 
SOC inside the triangle by assuming spin-flip hopping 
of strength a between each site, denoting also the SOC 
strength in analogy with the numerical model. Spin- 
uptown) corresponds to /i = 1 (— 1 ). 

The Hamiltonian of the system can be split in spin-up 
and spin-down blocks since the SOC does not mix the 


two spins, as 

H = 

where h f and /A are, 

( V 

= I 1 + i/ia 

y 1 — ifia 

by assuming real hopping t=l eV between the three sites, 
with /i = 1(—1) for spin up(down). The transmission 
probabilities between the different sites of the triangle 
can by calculated via the Green’s function 

G{E) = {EI-H-T,{E ))~ 1 (9) 

where E(E) = Ei (E) + T, 2 (E) + E 3 (.E) is the self-energy 
due to the semi-infinite linear chains at the corners de¬ 
scribed by H = JY c,q +1 + c.c with energy dispersion 
E = 2 tcos{k). E(EJ) is given by equal diagonal blocks for 
spin up and down 

/ E ld (E) 0 0) \ 

E "(E) = 0 Yii d {E) 0 (10) 

V 0 0 E ld (£) ) 

where T,i d {E) is the self-energy of each chain, which can 
be calculated recursively as follows [50]. 

In general we can calculate the Green’s function via 

G(E) = (E — g(E))~ 1 . (11) 

where g{E ) represents the surface Green’s function of the 
semi-infinite chain. Adding a single site via hopping t=l 
to the chain should not alter it’s Green’s function, so that 
G{E) = g{E). Consequently we can write the following 
recursive relation for g(E ), 

g(E) = (E — 3 (E ))- 1 =* g{E) = f ± i^l E, (12) 

Since we are interesting only in the retarded Green’s func¬ 
tion we choose the minus sign representing outcoming 
waves from the point of excitation, giving the self-energy 
E id{E) = g(E ), 

E I F2 

E ld (E) = - - i^J 1 - — for - 2<E<2. (13) 

Since H and E (E) are both block-diagonal in the spin 
basis, this is also true for Eq. [9] Therefore the Green’s 
function elements between opposite spins are zero, and 
we can write 

G»{E) = {EI -h*- E^(E ))- 1 (14) 

By plugging Eq. [H Eq. [TU] in Eq. [14] and inverting, 
we derive the following formula for the Green’s function 
element from site 1 (excitation) to 2 (response) 


/V 0 
0 


( 7 ) 



( 8 ) 




G^{E) 


1 -a 2 + E- Yii d (E) + ifia(- 2 + E - E ld {E) -V)-V 
(2 -E + E ld (E) + F)(-3a 2 + (-1 - E + E ld {E) + V ) 2 ) 


In order to calculate the respective transmission proba¬ 
bilities we apply the Fisclier-Lee relations, 

T&{E)=v{E?\G% 1 (E)\ 2 ( 15 ) 

where v(E) = ^ = —2^/1 — is the group velocity 
of the linear chains, assuming that h = 1. We note that 
the transmission probabilities between opposite spins are 
zero since the respective Green’s function elements are 
zero. 

By using the above relation we can calculate the spin 
transmission T( 2 (. E ) which is plotted in Fig. [G] 

T^(E)={tI 2 (E)~T^ 2 (E)). (16) 
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